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1. Introduction. 

The hadron production near heavy quark threshold will be thoroughly studied experimen- 
tally at future accelerators, e.g. [|I|]. The motion of the heavy quark- ant iquark pair near 
the production threshold is nonrelativistic to high accuracy that justifies the use of the 
nonrelativistic quantum mechanics as a proper theoretical framework for describing such a 
system [@, |^, Q. Being much simpler than the comprehensive relativistic treatment of the 
bound state problem with Bethe-Salpeter amplitude p, this approach allows one to take 
into account exactly such an essential feature of the near-threshold dynamics as Coulomb in- 
teraction 1^, 0] . For unstable heavy quarks with a large decay width it is possible to compute 
the cross section near threshold point-wise in energy because the large decay width serves 
as an infrared cutoff and suppresses the long distance effects of strong interaction [Q . 

2. tt-system near threshold 

The tt-pair near the production threshold is just a system that satisfies the requirement of 
being nonrelativistic. Therefore the description of tt-system near the production threshold 
y/s ~ 2?Tif {y/s is a total energy of the pair, nit is the top quark mass) is quite precise within 
NRQCD. Reasons for this accuracy are related to the large mass of the top: 

• The top quark is very heavy rrit = 175 GeV and there is an energy region of 
about 8 — 10 GeV near the threshold where the nonrelativistic approximation for 
the kinematics is very precise. The velocity is small (for E of interest in the range 
|E| ~ 2 ^ 5 GeV) 
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Relativistic effects are small and can be taken into account perturbatively in v (even 
in v"^). 

the strong coupling constant at the high energy scale is small asirrit) ~ 0.1 that 
makes the mapping of QCD onto the low energy effective theory (NRQCD), which is 



perturbative in asirrit), numerically precise. 

• The decay width of top quark is large, Tt = 1.43 GeV; infrared (small momenta) region 
is suppressed and PT calculation for the cross section near the threshold is reliable even 
point-wise in energy. 

The above properties could make the ti system near threshold flat and not interesting for 
physical study, however, there is a feature of quark-antiquark interaction that becomes domi- 
nant in this kinematical regime and brings a nontrivial structure into the dynamics. Because 
Q^s ~ f and the ratio Og/v is not small, the Coulomb interaction is enhanced. The ordinary 
perturbation theory for the cross section (with free quarks as the lowest order approxima- 
tion) breaks down and all terms of the order {ag/v)"' should be summed. The expansion for 
a generic observable f{E) in this kinematical region has the form 

f{E) = fo{as/v) + asfMs/v) + alh{a,lv) + ... (2) 

where fi{ois/v) are some (not polynomial) functions of the ratio Os/f , fo{ois/v) is a result 
of the pure Coulomb approximation (or a kind of its improvement). The expansion in a^ in 
eq. (0) takes into account the perturbative QCD corrections to the parameters of NRQCD 
and relativistic corrections (in the regime v ~ a^). One can study the ti system near 
threshold in the processes e+e~ -^ ti 0, [lO| and 77 — > ti [|rT| |l^. These processes have the 



following features: 

• e^e'^ — > ti: the production vertex is local (the electromagnetic current in case of 
photon and/or weak current in case of Z-boson), the basic observable is a production 
cross section which is saturated by S-wave (for the vector current), NNLO analysis is 
available. 

• 77 — > ti: the production vertex is nonlocal (T-product of two electromagnetic cur- 
rents), both S- and P-waves can be studied for different helicity photons, the number 



of observables is larger (cross sections as, crp, S-P interference). Because of nonlocality 
of the production vertex the high energy coefficient (necessary for mapping the QCD 
quantities to NRQCD ones) is more difficult to obtain. Some calculations were done in 
NLO [|l^, Q and the full analysis of the cross section is available in NLO of NRQCD 
only [|15]. The low energy part of the process can be studied in NNLO without a strict 



normalization to full QCD (see |T^ for relativistic corrections). 



3. Theoretical description. 

I shall discuss only e+e^ —>■ ti process mediated by the photon because NNLO analysis, 
which contains the most interesting features, is possible. The basic quantity is the vacuum 
polarization function 

n(^) = I j{Tj,^{x)jem{^)y'''dx, q^ = (2mt + Ef . (3) 

Near the threshold (for small energy E) NRQCD is used. The cross section is saturated with 
S-wave scattering. In this approximation the polarization function near the threshold to the 
NNLO accuracy in NRQCD is given by 

n(£;) = ^Ch{a,)Co{Elmt)G{E- 0, 0) . (4) 

mi 



The pole mass definition is used for rat (e.g. [0), «« is the strong coupling constant. 
The choice of normalization points for coupling constants entering different parts of the 
theoretical expression (^) will be discussed later. Ch[p^^ is the high energy coefficient which 
has been known in the NLO since long ago [IS, |19|, |2^ (before the explicit formulation 
of NRQCD). G(E;0,0) is the nonrelativistic QY , E = 4~s - 2mt. The leading term of 
the cross section in the effective theory representation is given by a correlator of currents 
with dimensionality 3 in mass units. The quantity Co{E /mt) describes the contributions of 
higher dimension operators within the effective theory approach. These contributions have, 
in general, a different structure than the leading term. However, to the NNLO of NRQCD 
the contribution of higher dimension operators can be written as a total factor Co{E /mt) 



for the leading order GF, Co{E/mt) = 1 — AE/3mt. The polarization function near the 
threshold (^) contains expansions in small parameters a^ and/or v, cf. eq. (0). The leading 
order approximation of the low energy part is the exact Coulomb solution for the Green's 
function. The more detailed description of the ingredients of the representation in eq. @ is: 

• the nonrelativistic Green's function G{E; x, x') is given by G{E) = {H — E)^^ where 

2 

H=^ + V(r) (5) 

mt 

is the nonrelativistic Hamiltonian describing dynamics of the tt-pair near the threshold. 

The most complicated part of Hamiltonian (P) to find is the heavy quark static potential 

Vpot{q) entering into the potential V{r). The static potential Vpot{q) is computed in 

perturbation theory 

%ot{q) = -Cf^ (l + asih In jj/q + ai) + a1{bl In^ /i/g + 62 In yu/g + 02) + . . .) • (6) 

Here 61 = 2/5o, 62 = 2(/3i + 2/5oai), 
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(4vr)/3o = -Ca - -Tpuf = - (7) 

is the first coefficient of the /9-function, 

(47r)2/5i = ^d - ^CATprif - ACpTpUf = H^ (8) 

is the second coefficient of the /3-function. The static potential can be written in the 
form 

VpoM) = -Cf^^ (9) 

that gives a definition of the effective charge ay related to the MS-scheme coupling 
constant 

«v(/w) = «s(Ai)(l + aiOsifi) + a2as{^if) ■ (10) 



Coefficients ai,2 are known. The numerical value for Oi reads [^ 
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The coefficient 02 is given by [^, |2^ 

,^ ,2 /4343 ^ 2 TT^ 22^,^, \ ^2 /1798 56^,^,\^^ 



V 162 

16C(3)) CfTfU} + (yT^Fn/) 
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V 3 
or in a more concise form, 

1 /7217 2 gvr^ 62C(3)\ 1 

02 = - — - + 367r2 ^^ = - — -155.842... 

(47r)2 \^ 162 4 3 y (47r)2 

In QCD we have Ca = 3, C^ = 4/3, Tp = 1/2, and n/ = 5 for the energy region of 
the tt-pair production. The effective couphng ay is nothing but a running couphng 
constant in some special subtraction scheme. The coefficient 02 allows one to find the 
effective /3- function jSy for the evolution of the coupling ay in the NNLO. Therefore 
the effective coupling constant for the static potential is now fully determined at the 
NNLO. 

High energy coefficient Ch{as) is given by the expression 

In NNLO, there appears a new term proportional to the logarithm of the factorization 
parameter fif that separates long and short distances (or large and small momenta) 
within the effective theory approach. The finite (/i/ independent) coefficient C2 is 



known 24, 25 
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The coefficient of NLO in eq. (|TT|) is /i/ independent - the factorization procedure 
(a separation of scales) is insensitive to the border. An explicit dependence of high 
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and low energy quantities on the factorization scale /i/ is a general feature of effective 
theories which are valid only for a given region of energy. A physical quantity, which 
is given by a proper combination of results obtained in different energy regions within 
respective effective theories, is factorization scale independent e.g. [|^, |^. In NRQCD 
this feature reveals itself in fif dependence of Green's function and of the high energy 
coefficient Ch- To see how this dependence emerges one can consider the contribution 
of a generic NNLO (three-loop) diagram into the cross section. In full QCD one can in 
principle compute it for an arbitrary numerical value of the invariant s = {2mt + E)"^ 
but it is a rather complicated function of the ratio s/rn^. It is difficult to obtain such a 



function analytically (numerical study was done in pH). However, the analytical result 
for the diagram in the threshold limit simplifies and leads to a logarithmic singularity 
in energy E (for the moment we forget about Coulomb singularities) 

a^ln(E/mt), E^O. (12) 

Within the effective theory approach this singularity splits between GF (the low energy 
singularity ln(£'//ij)) and the high energy coefficient Ch{ois) (the mass singularity 

ln(^/M)) 

al \Yi{E/mt) = al (HE/fif) + Hfif/nit)) . (13) 

In both limiting cases {E = 0, or calculation at threshold, and calculation within 
low-energy theory for the GF) the resulting integrals are simpler than the original 
integral for the diagram and can be done analytically. For one massive particle the 
calculation of diagrams is technically simpler and the complete expressions in full 
theory are available for some physical quantities ||2^. The corresponding effective 
theory near the static limit is HQET which is also simpler than NRQCD. In this case 
one can analyze the cancellation of factorization scale dependence between matching 



(high energy) coefficients ||3^ and the low energy (HQET) amplitudes explicitly. 



Collecting the above ingredients together allows one to have NNLO accuracy for Il{E) near 

the threshold. 

4. Solution for the low energy part. 

The main dynamical object used to describe the ti system near the threshold is the non- 

relativistic Green's function G = {H — E)^^. The Hamiltonian is represented in the form 



31, 32, 33 



2 

H = Hc + AH, Hc = -~ Cf— (14) 

m r 



with 



/\H = /\Vpot - 
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{^^'-t)m. (15) 
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Constructing the Green's function is straightforward and can be done analytically within 

perturbation theory near Coulomb Green's function Gc{E) or numerically (for complex 

values of E only that can be used to describe the production of particles with nonzero 

width) [^ |3^, ^, ^ |3^, |39|, ^ ^, 0. Results are presented basically as an expansion in 



consecutive orders 

G = Go + AGi + AGa (16) 

to check the convergence of the approximations. Here: 

• LO: Coulomb approximation, Gq = Gc 

• NLO: AGi -^ 0{as) corrections from the static potential Vpotij-) 

• NNLO: AG2 -^ 0{a1) corrections from a1 terms in the static potential Vpot{f^) and 
from the second iteration of the 0{as) term in Vpotir)., relativistic v'^ corrections. 

The relativistic H^ and anticommutator corrections in eq. ([T5|) can be taken into account 
by the shift of the parameters of the Coulomb Green's function E ^ E + E^/Arrit and 
tts — ^ as(l + 3E/2'mt). In this respect the modified Coulomb approximation can be used as 
the leading order approximation. The analytical solution for Green's function is perturbative 

G = Gc- GcAHGc + GcAHGcAHGc - • • • (17) 



The MS-scheme for the static potential Vpot{r) has been used in the solution. The numerical 



results obtained by the different authors are rather close to each other |53]. 
5. Features of the physical result. 

The top quark width Ft plays a crucial role in the calculation of the ti production cross 
section near the threshold ^. At the calculational level the width can be taken into account 
by a shift of the energy variable E. Operationally one can proceed as follows. The mass 
operator of the top quark is approximated by the expression M = rrit — iTt/2. Then the 
kinematical variable s — Amf relevant to the near-threshold dynamics is substituted with 
s — 4M^ {^/s = E + 2mt) and one finds 



s - 4M2 = Amt{E + iVt) + E^ + T 



t ■ 



Neglecting higher orders in E and Tt one obtains a recipe for taking into account the width 
r^ by the shift E —>■ E + iTf. The dispersion relation for the vacuum polarization function 
n(-E) has the form 

THE) - ' "f^''''^' 



E' ~E 
With the shift recipe one finds 

Because the point E + iTf lies sufficiently far from the positive semiaxis (and the origin) in 

the complex energy plane the cross section eq. (|18D is calculable point-wise in energy. For 

the bb system where F;, is small the situation is different and only moments of different kinds, 

for instance, 

p{s)ds 



are meaningful in the near-threshold Coulomb PT calculations. The hadronic cross section 
cr{E) was obtained by many authors (as a review see, |0). The normalized cross sections 
R'iE) for typical values m* = 175 GeV, F* = 1.43 GeV, as{Mz) = 0.118 are plotted in Fig.l 
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Fig. 1. The normalized cross section R'"[E) in the leading order (solid lines), NLO 

(bold dotted lines) and NNLO (bold solid lines) for rrit = 175 GeV, Tt = 1.43 GeV, 

as(Mz) = 0.118 and /t^ = 50 GeV, 75 GeV and 100 GeV. The dotted line corresponds 

to the result in Born approximation. 






^ . The curves have characteristic points which are usually considered as basic observables. 
They are: Ep - the position of the peak in the cross section and Hp - the height of the peak 
in the cross section. In the limit of the small Tt (at least for Tt which is smaller than the 
spacing between the first two resonances) one would have Ep ~ Eq and Hp ~ |'?/;o(0)p. 
For the actual value of Tt = 1.43 GeV which is larger than the spacing being the first two 
resonances the peak position and height are not determined by the first resonance only. The 
convergence for Ep and Hp in consecutive orders of perturbation theory near the Coulomb 
solution is not fast in the MS-scheme. For the typical numerical values of the theoretical 



parameters rrit, Tt and as{Mz) one finds [^ 



E„ 



Hr, 



Eo(l + 0.58 + 0.38 + ...) 
i7o(l- 0.15 + 0.12 + ...) 



(19) 



see also [^4|, 0). Important contributions that affect the quality of convergence are the 
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local term (~ asVoS{'r) which is related to 1/r^ non-Abelian term |^) and higher order PT 
corrections to Vpot{r). 
6. Possible improvements. 

The slow convergence for the peak characteristics of the cross section has caused some dis- 
cussion. The suggestions of the redefinition of the top quark mass have been made (as a 
review see [^ , also [^ , ^ ) . The use of the pole mass as a theoretical parameter for the 



description of the cross section near the threshold is criticized on the ground of its infrared 
instability [^. It is usual that some effects of interaction can be taken into account by 



introducing the effective mass parameter for the particle 1 50, 51, 52]. In this talk I only 
discuss some possible ways of optimizing the convergence for the Green's function with the 
pole mass as a theoretical parameter. Note that actual calculations near the threshold have 
been performed within the pole mass scheme. For optimizing the convergence one can use 
methods of exact summation of some contributions in all orders and renormalization scheme 



invariance of PT series e.g. |3J]. The Hamiltonian can be written in the form 



H = HLo + AVpT + asVoS{f), Vo = -^(^ + ^) (20) 

where corrections are given by the perturbation theory corrections to Vpot{r) (AVpy-part) and 
by the local term (asVo(5(r)-part). These two contributions can be dealt with more accurately 
than in the straightforward approach. Indeed, the (5(f^-part is a separable potential and can 
be taken into account exactly |5^. The solution reads 



with 

GUE) = {Hlo + ^VpT - E)-' (22) 

being the irreducible Green's function. Dealing with the PT expansion of the static potential 
in NRQCD is an important issue in getting stable results for the cross section near the 
threshold because the static potential is the genuine quantity which is computed in high 
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order of PT in the strong coupling constant ^EM . The convergence in the MS scheme is not 



fast which reflects the physical situation that the observables represented by the cross section 
curve (for instance, Ep and Hp) are sensitive to different scales. The finite-order perturbation 
theory expansion of the static potential given in eq. @ cannot handle several distinct scales 
with the same accuracy. Indeed, the PT expansion of the static potential is done near some 
(arbitrary) scale (or distance) which can be considered simply as a normalization point. The 
farther a given point lies from this normalization point the worse the precision of the PT 
expansion for the static potential at this point is. The PT expressions in the MS scheme are 
not directly sensitive to physical scales because subtractions are made in a mass independent 
way (for instance, massive particles do not decouple automatically in the MS scheme e.g. 
||55|). Therefore it is instructive to rewrite the static potential in more physical terms than 
just the MS-scheme parameters (aj^, /i) (which would remind the reader of the momentum 
subtraction scheme where, for instance, the decoupling is explicit [^) 



Vpotir) = -C^— (l + «o&ilnr/ro + a^(6?lnV/ro + 6'2lnr/ro + c) + ...) . (23) 



Here r^ and c are the parameters of the renormalization scheme freedom in NNLO and oq is 
the corresponding coupling in the {ro, c}-scheme [^. They parameterize the center of the 
expansion (a normalization point) and the derivative (respective /3-function) of the static 
potential. The parameters (ro, c) can be chosen such in order to minimize the higher order 



corrections to a particular observable (e.g. [Q where NLO analysis has been done). In such 



a case ro can be understood as a typical distance to which a chosen observable is sensitive. 
Note that the best approximation of the static potential Vpot{r) for different scales would be 
provided by the use of the running coupling constant as{r). The analytical calculation of 
the Green's function becomes just impossible in this case. However, even for the numerical 
calculation of the Green's function one cannot naively use the running coupling constant 
as(r) in the static potential Vpot{r) for all r. With a generic running coupling constant, 
the IR singularity (Landau pole) can occur in as{r) at large r which are formally necessary 
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for computing the Green's function. This singularity has httle to do with usual divergences 
within the effective theory approach and cannot be removed by the standard renormalization 
tools. It can be dealt with if an IR fixed point appears in the evolution for the effective 



coupling constant (e.g. [^) or with some other regularization e.g. |^9[. For the top quark 
production the contribution of this area (large r) into the cross section is small because of 
the large decay width of the top quark. In the finite-order PT analysis the parameters ro and 
c can be chosen to minimize higher order corrections either to Ep or to Hp but not to both 
simultaneously because Ep and Hp are sensitive to different distances. Indeed, one finds the 
difference of scales minimizing corrections to the first Coulomb resonance in NLO to be 

Hte/t^) = \ + \- (24) 

Because of the large top quark width many states (resonances and continuum alike) con- 
tribute into the position and height of the peak in the cross section. Therefore the charac- 
teristic distance estimates are not so transparent (the NNLO peak position, for instance, is 



not exactly the ground state energy in zero width limit). The relation (p^) can serve just as 

a basic guide. In practical analysis one can choose the particular numerical values for the 

parameters (rg, c) which stabilize either Ep or Hp. 

7. Conclusion 

To conclude, we have presented a cross section for the ti pair production near the threshold. 

The result is based on the solution to the Schrodinger equation for the nonrelativistic Green's 

function. Some methods of resummation and convergence optimization are discussed. 
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